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i5 ■ ABSTRACT 

£\j | We recently proposed that molecular cloud dense cores undergo a prolonged period 

of quasi-static contraction prior to true collapse. This theory could explain the observa- 
tion that many starless cores exhibit, through their spectral line profiles, signs of inward 
motion. We now use our model, together with a publicly available radiative transfer 
code, to determine the emission from three commonly used species - N2H + , CS, and 
HCN. A representative dense core of 3 M that has been contracting for 1 Myr has line 
profiles that qualitatively match the observed ones. In particular, optically thick lines 



Oh- 



o 

have about the right degree of blue-red asymmetry, the empirical hallmark of contrac- 
tion. The J = 2 — >• 1 rotational transition of CS only attains the correct type of profile 
if the species is centrally depleted, as has been suggested by previous studies. These 
results support the idea that a slow, but accelerating, contraction leads to protostellar 
collapse. In the future, the kind of analysis presented here can be used to assign ages 
\Q . to individual starless cores. 

Subject headings: ISM: clouds, kinematics and dynamics, molecules — stars: formation 
Q \ — line: profiles 



^ ■ 1. Introduction 
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Over the last decade, astronomers have carefully scrutinized the large fraction of molecular 
cloud dense cores that contain no point sources of radiation, in order to understand more fully 
the onset of star formation. It was already established t hat the gross prope rties of these objects 



differ little from their counterparts with i nternal stars (|Lee fe Myerd I1999T ). except perhaps for 
a mass that is slightly higher on average ( Jijina. Myers. Sz Adamslll999l . Fig. 3). More detailed 



spectroscopic studies of molecular line em i ssion revealed that many s t arless cores exhibit signs 



of inward contraction ([Williams et al.lll999l : iLee. Myers. &: Tafallalll999l : iGregersen fc Evan 
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Lee. Myers. &: Tafallal l200ll ; ISchnee et al.l 120071 ) . The hallmark of this so-called infall signature is 
an asymmetric emission line, often with a self-absorption dip, that is skewed to the blue. 

There is a variety of line profiles seen in dense cores, sometimes even in the same object, 
and the profile just described is by no means universal. It is pre ferentially found in starless 
cores with a higher cent ral column density and optical depth (e.g. iLee. Myers. &: Tafallal 1 19991 : 
Gregersen fc Evansl [2000) . Such cores also exhibit chemical abundances that differentiate them 



from lower-density objects (jCrapsi et al.ll2005l : iKeto fc Casellil l2008l ) . The kinematic signature is 
especially common in m ore evolved dense cores that co ntain Class and Class I sources, i.e., very 
young, embedded stars (JMvers. Evans, k, Ohashill2000l ). These facts, togeth er with the simple and 
compelling infall interpretation of the asymmetric profiles (JMvers et al.lll996l ). have convinced most 
researchers that we are witnessing a type of contraction that precedes the free-fall collapse onto a 
protostar. The contraction hypothesis is especially convincing for starless cores, where there is no 
possibility of confusion from stellar outflows. 

What is the physical cause of this motion? The observed magnitude of the infall velocity is 
typically a few tenths of the sound speed. This finding rules out ambipolar diffusio n as the under- 
lying process, since the characteristic ion-neutral drift speed is much smaller (e.g. ICiolek fc Basu 



2001 



) . Conversely, t he subsonic level of the velocity, together with the statistical prevalence 



of starless cores (e.g. 



Forbrich et al. 



20091 ). indicate that such objects are relatively long-lived 



( Kirk. Johnstone, fc D i Fran cesco! 20051 ). and thus are unlikely to be in a true state of collapse, 



as proposed by lKeto fc Field! (120051 ). 



In a recent contribution (jStahler fe Yenll2009l . hereafter Paper I) we introduced a new model for 
the contraction. Starless cores, like all self-gravitating objects, undergo bulk oscillations. Indeed, 
there are now several cases where spectros c opic observations ind i cate that such oscilla tions are 
occurring (IRedman. Keto. fc Rawlingd 120061 : iBroderick et al.l 120071 : IVelusamy et al.l 12008) . For an 
object that has evolved to the point of collapse, by whatever means, the frequency of the lowest 
normal mode vanishes. We studied, using perturbation theory, how a cloud subject to such a frozen 
(more commonly called 'neutral') mode slowly evolves toward full-blown collapse. We found that 
this previously unrecognized phase of quasi-static contraction lasts a considerable period (of order 
10 6 yr), during which time the cloud attains internal speeds similar to those observed (see Fig. 4 
of Paper I). 

Encouraged by this result, we now seek to forge a stronger link to the observations. We 
calculate, using a publicly available radiative transfer code, the emitted line spectrum from our 
model cloud. We choose three molecules that are among the most commo nly employed in the 
spectr osco pic studies: N2H" 1 ", CS, and HCN. Following the early surveys of iTurner fc Thaddeus 
(J1977I ) and lWomack. Ziurvs. fe Wvckoffl (J1992I ). N2H + is often used as a tracer of dense, relatively 
quiescent cloud gas. Its rotational lines are optically thin in low-mass dense cores. Thus, they do not 
exhibit either self-absorption or blue-red asymmetry, but serve principally as a gauge of excitation 
temperature and as a comparison benchmark with lines that are sensitive to internal cloud motion. 



The rotational lines of CS are ofte n optically thick, an d can have asymmetric profiles that have 
long been interpreted kinematically (jWalker et al.lll986h . On the other hand, the infall signature is 
relatively weak in starless cores, in par t because the polar mo lecule CS is depleted at high density 
from the gas phase onto grain mantles (JBergin fc Langerlll997l ). Finally, the J = 1 — y rotational 
transition of HCN has the strongest infall signa ture of all, and the molecule appears to suffer little 



depletion (JAikawa et al. 



2005 



Sohn et al.ll2007l ). 



Although our model cloud is a highly idealized isothermal sphere, the results of this study 
further lend credence to the underlying dynamical picture. For a representative 3 M Q cloud that 
has been contracting for 1 Myr, we find that all the calculated line profiles resemble those observed, 
at least in an average sense. Interestingly, the CS profile is only acceptable once we include central 
depletion of the molecule. In Section 2 below, we first describe our implementation of the radiative 
transfer program. Section 3 then summarizes observations to date of the three molecules in starless 
dense cores. We describe, in Section 4, how the emergent line profiles alter when we change our 
physical input parameters. Judging from these trends, we then select a canonical model cloud and 
compare its emission spectrum to the available data. Finally, Section 5 discusses limitations of the 
dynamical theory and the prospects for improvement. We also indicate how the model, even in its 
current, simplified form, can be used to assign contraction ages to individual starless cores. 



2. Calculating Line Emission 

For ea ch selected cloud model, we solve the equation of transfer numerically, using the code 
RATRAN (jHogerheijde fe van der Tak!l2000l ). The basic strategy underlying this code is to choose 
randomly a set of cells outside the cloud that generate photons. These photons propagate through 
the cloud interior, while the molecular level populations are temporarily held fixed. The statistical 
equilibrium of the levels is then solved separately, whe reupon photon propagation is reinitiated. 
The accelerated lambda iteration scheme introduced by lRybicki &: Hummed (|199ll ) ensures rapid 
convergence of the radiation field and level populations. Treatment of non-isotropic scattering is 
cumbersome in this approach, so RATRAN does not include this effect. Such a simplification is 
acceptable in the infrared and millimeter regime of interest. 

The code must be supplemented by tables of molecular data, containing energy levels, degen- 
eracies, Einstein A-values, and c ollisional rate coeffi cients. Here, we utilized the Leiden Atomic and 
Molecular Database (LAMDA). IShoier et al.l (120051 ) describe the construction of the database and 
its use. These tables include the full multiplets of rotational lines created by hyperfine splitting. 
In fact, recent observational studies of dense core kinematics have focused on individual hyperfine 
lines, as we will describe. 



In order to display the line spectra, we utilized the MIRIAD software package (jSault. Teuben. &: Wright 



19951 ). We placed our cloud at a representative distance of 140 pc, and "observed" it with a Gaus- 
sian beam of angular size 60". For comparison, our canonical cloud has a physical diameter of 



0.27 pc, or 392" at the selected distance. In constructing spectra, we assumed a velocity resolution 
of 0.01 km s -1 . 

We tested the radiative transfer calculation by first repro ducing prior results for line profiles 
in a collapsing cloud. An early and influential work was that of lZhoul (119921 ). This s tudy compared 
the self-similar col la pse models of IShul (|1977l ) with that derived independently by iLarsonl (|1969) 
and iPenstonl (|1969l ). IZhoul (|1992D calculated, for both models, CS rotational line profiles at various 
times, both before and after formation of the central protostar|^j 

In testing bot h self-similar collaps e models, IZhoul ( 19921 ) assumed a gas kinetic temperature 
of 10 K. Following iLizano fe Shul (|1989l ). he augmented the Shu model with a turbulent velocity of 
the form 

Vturb = \[kTp , (1) 



where n is a constant and p is the local cloud density. The constant was d eterm i ned em pirically by 
Myers fc Goodman! (j 19881 ) to have the value k = 1.1 X 10~ n g cur 1 s~ 2 . IZhoul (|1992l ) argued that 
this figure was likely to be an overestimate, since the observed "turbulence" could have included 
infall motion for dense cores containing stars. Accordingly, he used a value of k smaller by a factor 
of 0.25. Finally, he set the CS abundance, i.e. , the number density of the molecule relative to 
hydrogen, to 4 x 10~ 9 , following iFullerl (|1989l ). 



We provisionally adopted the same assumptions as IZhoul (|1992l ) and recalculated, using RA- 
TRAN, CS line profiles from the Shu model for the period after protostar formation. All lines were 
viewed through the cloud center. The RATRAN code allows one to select N tot , the total number 
of shells in the cloud, as well as r out /ri n , the ratio of outer to inner shell radii. We adjusted both 
numbers until we obtaine d the best m atch with the J = 1— > 0, J = 2 — > 1, and J = 3 — > 2 pro- 
files shown in Figure 3d of lZhoul (J1992I ). By choosing N = 100 and r„ ^t /r\ n = 500, we obtained an 
excellent fit to all profiles for the four times considered by lZhoul ([1,992 ). Agreement did not improve 
when either N tot or r out /ri n was increased further. For example, the peak antenna temperatures 
Ta changed by less than 0.1 %. We therefore retained these values of A^tot and r ont /r- m in our 
subsequent calculations. 

Turning to starless dense cores, we continued to use a gas kinetic temperature of 10 K, corre- 
sponding to a therm al speed of ar = 0-2 kms" 1 . T his temperature is on the low end of the range 
found empirically by lJijina. Myers. fe Adamsl (|1999l . Figure 2). In any case, it is important to note 
that purely thermal line pro files are never seen, even at the centers of the most quiescent objects 
(JBarranco fc Goodmanlll9 98). and the prescription for turbulence in equatio n (1) continu es to pro- 
vide a useful representation of the observations. In light of the comments by IZhoul (| 19921 ) . we view 



1 The original Larson-P enston model only considered collapse up to the point of proto star format ion. IZhoul 1 1992 ) 
utilized the calculation of lHunterl (|l977T I to extend the flow past thi s point. Conversely. IShul (|l977T l only considered 
times after protostar formation. For the earlier epoch, IZhoul ( 1992 ) employed a sequence of hydrostatic spheres of 
increasing central density. 



k as a free parameter, whose upper limit is likely to be the value advocated bv lMvers &: Goodman 
(|1988l ). It is straighforward to incorporate turbulence into RATRAN, which allows a variable degree 
of intrinsic line broadening in each mass shell of the cloud. 

An additional free parameter is inherent in the theoretical model of Paper I. We envisioned 
the cloud undergoing a global oscillation, characterized by a nondimensional amplitude e. For our 
standard initial conditions, e gives the fractional deviation of the density from its unperturbed, 
equilibrium distribution (see eq. (17) and the discussion in §4.2). The numerical integration of 
Paper I used a re latively large am plitude of e = 0.2, motivated by observation and analysis of the 



starless core B68 (Ketoetal 



20061 ). Here, we reset e to the more modest, and perhaps more typical, 



value of 0.1. We comment later on how changing e affects the emergent line profiles. 

With these assumptions, the only remaining free parameters are the cloud mass M and the time 
of observation t. The zero point for t is the onset of quasi-static contraction. Both quantities must 
be related to the corresponding nondimensional ones, M and t, used in Paper I. Our model cloud 
has a fixed nondimensional mass, M = 4.19, where the number results from numerical integration 
of the isothermal Lane-Emden equation. Combining equation (9) and (12) of Paper I, we find that 
M and t are related by 



Ma\t 
GM 



(2) 



Since M and a? are either known or preset, the equation tells us, for any M and t, the appropriate t 
to use from the model sequence. We further read off, at this point in the sequence, the distribution 
of the nondimensional density and velocity as functions of radius. The equations of §2.2 in Paper I 
tell us how to translate the quantities into dimensional formo 



Observations of the Molecules 



3.1. N 2 H+ 

After surveys had established the presence of N2H + in star-forming r egions, targeted observa- 
tions were pursued in dense cores, inc l uding those already detected in NH3 (jCaselli. Myers. fe Thaddeus 
19951 ; Benson. Caselli. fe Mverdll998l ; ICaselli et al.ll2002l ) . It was found that the molecule selectively 
traces dense and quiescent gas, and is indeed absent in hotter and shocked material. The J = 1 — > 
rotational transition, which has 7 hyperfine components, is usually optically thin in starless cores, 
but not always so. For example, the well-studied core L1544 sho ws an asymmetric p r ofile toward the 



cloud center, extending over a region some 1400 AU in radius (JTafalla et al. 
1999n . 



1998 



Williams et al. 



2 These equations contain the external pressure P e . Equation (9) of Paper I allows us to solve for P E , given the 
dimensional cloud mass M. 



Lee. Myers, fc Tafallal (|1999l ) conducted the first largescale, systematic observations of N2H H 



in starless dense cores. The 220 sources in their catalog had been selected primarily on the basis 
of optical obscuration. For the 72 objects detected in N2H + , detailed study was made of the 
relatively isolated F\,F = 0, 1 — >• 1,2 hyperfine line at 93.176 GHz. This line almost always has 
a Gaussian profile, with no self-absorption or asymmetry. The mean measu red peak intensity was 



T A = 0.19 ± 0.11 K, and the FWHM was AV N2H + = 0.35 ± 0.14 km a -1 . iLee. Mvers. fc Tafalla 



(J200ll ) followed up these single-pointing observations with spatial maps of many individual cores. 
The N2H + emission is compact, with a well-defined central peak. Diameters of the regions within 
the half-maximum contours range from 0.05 to 0.13 pc. 



3.2. CS 



These properties of N2 H" 1 " emission are in sharp con trast to those of the J = 2 — ► 1 rotational 
line of CS at 97.981 GHz. ILee. Mvers. fc Tafallal (]1999l ) found, in the 163 starless cores for which 
the line was detected, that a large fraction of the profiles were asymmetric and skewed toward 
the blue, either with two distinct peaks or, more commonly, a blue peak and a red shoulder (see 
their Figure 1). Detailed profiles vary from source to source, and there are even a few examples of 



doubl e-peaked profiles stronger in the red. The subsequent mapping survey of lLee. Myers, fc Tafalla 



(J200ll ) showed that the region of CS asymmetry is relati vely diffuse, and the emission is frequently 
depressed just where the N2H + intensity peaks (see also lTafalla et al.ll2002l ). 



With both the optically thin (N2H + ) and optic ally thick (CS) l i nes in hand, one may quan- 
tify the degree of blue-red asymmetry of the latter. iMardones et al.l (j 19971 ) introduced SVqs, the 
normalized velocity shift of CS. This quantity is defined as 

5V CS = VcS A ~ ^ H+ . (3) 

AI/ N2 h+ 

Here, Vqs and Vn 2 h+ are the peak, line-of-sight velocities of the two species. The peak velocity of 
N2H + is assumed to represent the syst ematic motion of the cloud as a whole. In the single-pointing 
survey of lLee. Mvers. &: Tafallal (J1999I ). SVqs was usually negative, signifying that the profile peak 
is blueward of line center. The mean value of SVcs hi this sample was —0.24 ± 0.04. 



3.3. HCN 



The J = 1 — > rotational transition of HCN is the latest diagnostic tool for studying the 
kinematics of starless cores. As mentioned previously, the mo lecule suffers little depletion, at 
least for central densities less than 10 5 cm -3 ( Aikawa et alj|2005l ). Integrat ed intensity map s show 
that the HCN emission region is roughly coextensive with that of N2H + (jSohn et al.ll2004l ). The 
J = 1 — > transition consists of three hyperfine lines of differing optical depth: F (0 — 1) at 
88.634 GHz, F(2-l) at 88.632 GHz, and F(l-l) at 88.630 GHz, Under optically thin, LTE 



conditions, these would have intensity ratios of 1 : 5 : 3. Of course, it is the fact that these lines 
are optically thick in starless cores that makes them useful kinematic probes. 



Sohn et al.l (|2007l ) conducted single-pointing observations for HCN in 85 starless cores, netting 



64 detections. The profiles of all three hyperfine lines are distinctly non-Gaussian, displaying either 
two peaks with a self-aborption dip or a single peak with a shoulder. Not surprisingly, the central 
dip is strongest in the optically thickest F (2 — 1) line, which also shows the greatest blue-red 
asymmetry. But a central result of this study is the heterogeneity of the profiles. While 43 percent 
of the sources have at least one or two double-peaked hyperfine lines skewed to the blue, 33 percent 
have blueward and red ward asymmetries mixed in the three lines. 

Nevertheless, the overall tende ncy is for pr o files i ndicative of contraction, rather than expan- 
sion. In analogy with equation (3), ISohn et al.l (J2007I ) defined 5Vhcn ; a normalized velocity shift 
with respect to N2H + . For each hyperfine line, the large majority of sources are observed to have 
negative SVucn (see their Figure 4). The actual magnitudes of 5Vucn are larger than for CS. 
The mean values are -0.51 ± 0.06, -0.41 ± 0.10, and -0.59 ± 0.08 for F(0-1), F(2-l), and 
.F(l — 1), respectively. Moreover, the magnitudes are tightly correlated, in the sense that a source 
with greater |<5Vhcn| i n the F (0 — 1) line also has a greater shift in F (2 — 1) and F (1 — 1). Within 
each source, however, there is no clear pattern for the relative magnitudes of the velocity shifts for 
the three lines, nor for the relative intensities of the lines. In other words, starless cores frequently 
exhibit HCN "intensity anomalies" with res pect to LTE, a phenomenon that has long been noted 
and studied in other molecular clouds (e.g.. iGonzales- Alfonso fc Cernicharolll993l ). 



Because of the range in optical depths in the three lines , the profi l es of HCN are potentially 
useful for tracing the full velocity structure of a dense core. iLee et al.l (J2007I ) constructed single- 
dish spatial maps of emission for two starless cores, L694-2 and L1197. To complement their 
observations, they also performed a radiative transfer calculation in spherical cloud models with 
adjustable density and velocity pro files. Inter e stingl y, monotonic velocity profiles do not match 
the data. For both starless cores, ILee et al.1 (|2007l ) found that the inferred contraction speed 
begins relatively small at the cloud edge, reaches a maximum about halfway inward in radius, and 
thereafter declines toward the cloud center. This type of velocity profile emerges naturally from 
our dynamical model (Paper I; Figure 4). 



Numerical Results 



4.1. Trends in Calculated Line Profiles 



After setting the oscillation amplitude to e = 0.1, we next gauge how the line profiles, still 
viewed through the cloud center, change as the other parameters are varied. These trends enable 
us to seek a model whose lines best match the data in an average sense. Recall that our remaining 
free parameters are: n, the turbulence coefficient in equation (1), the cloud mass M, and the time t 



since the start of inward motio n. The updated molecular abun dances we used in the transfer code 
were : X (N 2 H+) = 7 x KT 10 teenson. Casern fc Mverslll998h: X (CS) = 4 x 10" 9 (JTafalla et al 



20021 ); and X (HCN) = 3 x 10" 9 ((Fuller et al 



1991 



Li et al. 



2004) 



We first explored the effect of altering the cloud mass at fixed k and t. Following lLee. Myers. &: Tafalla 
(jl999l ). we focused on the F x , F = 0, 1 ->• 1, 2 hyperfine line of N 2 H+. i As M is raised, the N 2 H+ 



line remains Gaussian, and is therefore still optically thin. Of course, the total emission increases 
with M, so that AVn 2 h+ a ^ so rises. 

In CS and HCN lines exhibiting self-absorption, these dips eventually disappear with increasing 
cloud mass. This tendency is easy enough to understand. According to equation (2), an increase 
of M at fixed t necessitates a lower t. We are viewing our fiducial, contracting cloud at an earlier 
nondimensional time, when the cloud radius f was larger. But equations (9), (10), and (12) of 
Paper I imply that 

_ 2 rht a,T 



MR =—7T + 



(4) 



Since i/r 2 falls, the column density through the cloud decreases with increasing M. Many observed 
lines, especially those of HCN, show clear self-absorption. Thus, the sensitivity of the calculated 
profiles effectively sets a ceiling on the allowable range of cloud masses. 

A similar disappearance of the self-absorption dips in CS and HCN occurs when we fix M and 
t, but increase k. A larger value of k corresponds to a greater velocity dispersion throughout the 
cloud. The gas consequently has a higher, overall excitation temperature, so there is less of the 
absorption that creates central line dips. By the same token, the optically thin N 2 H + line broadens, 
i.e., AF N2H + increases. Both trends help us to set bounds on k. 

Finally, we may fix M and k, and consider various evolutionary times t. At greater t, the 
maximum infall speed within the cloud increases. Thus, the blue/red asymmetry of all optically 
thick lines becomes more pronounced. Quantitatively, both SVcs and SVucn which are both 
negative, become larger in magnitude. Figure 1 shows that the overall emission in N 2 H + increases 
with time, a reflection of the higher cloud density. The line peak grows faster than its width, with 
the result that <5Vn 2 h+ declines. In the figure, the FWHM is 0.24 km s _1 at the initial instant, and 
0.21 km s _1 at the last time shown (t = 1 Myr). 



4.2. Canonical Model 

Rather than attempt to find, within the model parameter space, matches to individual dense 
cores, we sought first a combination of M, t, and k that reproduces the general, observed charac- 



3 The N 2 H+ table of LAMDA lists 3 lines at 93.176 GHz corresponding to the transitions Fi,F = 0, 1 -> 1,2, 
Fi,F = 0, 1 — > 1, 1 , and F\,F — 0, 1 — > 1,0. The lower levels are also listed as being degenerate in Table 1 of 



L, ■ 



'' 



Daniel et al.l (J2005). Assuming the three individual lines cannot be resolved in practice, we added their intensities. 



teristics of the various line profiles, as summarized in §3. Some experimentation was required, in 
light of the constraints just given. However, a suitable parameter set quickly emerged. If we set 
M = 3 Afo, t = 1 Myr, and k = 5.5 x 10 



12 g cm s 2 (i.e., 0.5 times the original value found by 



Myers &: Goodman! (J1988I )). then the profiles are generally consistent with the data. 



Following our previous discussion, it is the magnitude of the self-absorption dip that is most 
sensitive to variations in these parameters. We may accordingly gauge the acceptable range of 
parameter values by noting when the dips depart substantially from those in the canonical model. 
Holding t and k fixed while allowing M to range from 2.7 to 3.2 M@, the magnitude of the dip in 
the F (2 — 1) line of HCN changed by less than 25 % in either direction. By the analogous measure, 
t can vary from 0.8 to 1.2 Myr, and k from 3.3 to 8.8 x 10 -12 gems -1 . 

The topmost curve in Figure 1 shows the calculated N2H + profile for the canonical model. As 
in nearly all observed cases, the line is symmetric and Guassian. Our FWHM of 0.21 km s _1 is on 
the low end of the observed range, as would be expected given our assumed kinetic temperature 
of 10 K. Conversely, our TA-value of 0.33 K is higher than the average. Here, of course, the beam 
efficiency (assumed to be unity in RATRAN) comes into play, as do the various distances to the 
observed dense cores. 

Figure 2 (dashed curve) shows the emergent profile of CS through the center of the canonical 
cloud model. Because of higher optical depth, the line is markedly asymmetric, as the data also 
show. However, while the sense of the asymmetry in t his calculated profile is correc t, the deep, 
self-absorption dip is rarely seen (compare Figure 1 of iLee. Myers. <fe Tafallal (|1999l )). What is 
observed instead is a relatively broad, redshifted shoulder joining onto a blueshifted peak. 

We may recover this shap e by introducing dep letion of the molecule, presumably caused by 
adsorption onto grain surfaces. iTafalla et al.1 ( 2002 ) invoked such depletion to improve the obser- 
vational match of their CS spatial intensity maps. Following these authors, we assume a centrally 
depressed, radial profile for the relative number density: 



X(r) 



X exp[-n(r)/n d ] 



(5) 



Here, X is the usual interstellar value of X (CS), n(r) is the cloud number density, and the threshold 
density n^ is a free par ameter. For t h e soli d profile in Figure 2, we used n<j = 1 x 10 cm -3 . As a 
comparison, Table 4 of ITafalla et al.l ( 20021 ) lists n^- values ranging from 1 to 5 x 10 4 cm -3 for five 
starless cores. Our calculate d relative velocity shift, <5Vcs = —0.24, is equal to the mean observed 
(|Lee. Mvers. fc Tafallal Kwl ). 



The three hyperfine lines of HCN (Figure 3) exhibit more diverse profile shapes. The optically 
thinnest F (0 — 1) transition has a blue peak/red shoulder profile, similar to that in CS. The 
other two lines exhibit clear self-absorptio n, with the deepes t trough being found in the optically 
thickest, F (2 — 1) line. Scanning Figure 1 of lSohn et al.l (120071 ). we see that this pattern is commonly 
seen (e.g., L63, L1197), but is by no means universal (e.g., L1544). Our calculated relative velocity 
shifts, 5V, for the (0-1), (2-1), and (1-1) lines are -0.33, -0.81, and -0.67, respectively. As mentioned 



10 



earlier, observations give similar values for all three lines, with a mean of about -0.5. We are further 
encouraged by this qualitiative agreement. 

We stress that calculated line profiles are quite sensitive to the interior dynamics of the cloud, 
and thus do provide a worthwhile test of the underlyi ng; theory. To further illustrate this point, 



we recalculated all profiles using the collapse model of IShul (j 19771 ). We again set t = 1 x 10 6 yr, 
where the zero of time now refers to the start of inside-out collapse. For T = 10 K, the mass of the 
central protostar in the model is 2 M & . In calculating profiles, we used the same velocity dispersion 
parameter, k = 5.5 x 10 -12 g cm s -2 , as in our canonical model. 

All calculated line profiles are now far broader than those observed. For example, the N2H 4 " 
profile, while remaining symmetric, is distinctly non-Gaussian, with a shallow, central depression 
and wings extending to ±lkm s _1 on either side of line center. Optically thick lines have profiles 
that are much too asymmetric. Our point here is not that the Shu model is a poor representation of 
starless cores, which was to be expected. The message is rather that the gas motions induced by a 
central protostar of solar-type mass are qualitatively different from those in our slowly contracting, 
starless cloud, and in a manner clearly reflected in the molecular line profiles. 



4.3. Determining e and t 

Our canonical value of the oscillation amplitude e, while guided by the scant observational 
data available, is still rather arbritrary. How, in the future, can we pin down this parameter when 
we examine specific dense cores? Even more importantly, will we be able to gauge the cloud's 
contraction age through comparison between real and model-generated line profiles? 

These questions are related. We remarked in Paper I that an increase of e, if accompanied by 
a decrease in the evolutionary time, could leave the cloud largely unchanged. If that were the case, 
it would be difficult to ascertain e and t individually. Specifically, our cloud model, at a given e and 
nondimensional time t, attains nearly the same contraction velocity as one at a different amplitude 
e' and nondimensional time t! = t (e/V) . The nondimensional time scales as the inverse square of 
the amplitude because cloud contraction is driven by a second-order perturbation. The dimensional 
time may be obtained from the nondimensional one through equation (2). Does this scaling relation 
mean that line profiles are degenerate in e and tl Armed with RATRAN, we now demonstrate that 
this is fortunately not the case. 

Figure 4 shows in detail the effect of trading amplitude and age. The solid curve in the lefthand 
panel is the velocity profile in our canonical model (e = 0.1; t = 1 Myr). The dashed curve is the 
profile for a cloud with the same mass of 3M , but with e = 0.2 and t = 0.25 Myr. As predicted, 
the velocities are quite similar. On the other hand, the internal densities differ greatly. The 
righthand panel shows that the cloud with larger e has a lower central density, since it was initially 
more inflated and evolved over a shorter time. 
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The net result is that the line profiles also do not match. As an example, Figure 5 displays the 
three hyperfine lines of HCN. The degree of blue/red asymmetry in all the profiles is markedly less 
than in the canonical model. For the optically thinnest, F (0 — 1) line, the profile is nearly Gaussian. 
The normalized velocity shift 5V is only —0.16 km s _1 , well below typical observed values. 

What happens if we keep e = 0.2, but increase t back to 1 Myr? For our model to be useful 
observationally, the resulting line profiles should again differ noticeably from the canonical case. 
Indeed they do, and in the opposite sense as the previous example. Both the CS and HCN profiles 
now exhibit too much asymmetry. For example, the red wing of the F (0—1 ) line of HCN degenerates 



to a s houlder; a similar profile is seen in CS. According to Figure 1 of iMvers. Evans. &: Ohashi 



(|2000l ). this development reflects relatively high infall speed. In fact, the maximum infall speed for 
this choice of parameters is 0.06 km s _1 , three times that in the canonical one. The central density 
is also higher by about the same factor. Consequently, the N2H + line, while maintaining a Gaussian 



profile , has a T^-value of 0.85 K, significantly above the range reported bv lLee. Myers. fe Tafalla 



tea ). 



4.4. Abundance Gradients 

We have seen that the calculated CS profile is not even qualitatively correct unless we include 
the effect of depletion onto grain surfaces. Our strategy, following previous authors, was simply to 
introduce a spatial variation in the molecular abundance. Nevertheless, our canonical model still 
assumes a uniform distribution of N2H + and HCN. This simplification also bears examination. 

The detailed abundance of any species is influenced not only by the efficiency of grain ab- 
sorption but also by the clo u d dyn amics itself. The latter point was emphasized in the pioneering 



study by iRawlings fc Yates! (|200ll ). who followed a time-dependent chemical reaction network in 
models of collapsing dense cores. Utilizing a radiative transfer code, these authors demonstrated, 
as we did in §4.2, the sensitivity of emergent molecular line profiles to the cloud's internal velocity 
structure. Not only is there a direct, kinematic effect, but various species may not have time to 
achieve chemical e quilibrium, creating abundance gradien ts that also influence line profiles. In a 



more recent work, iTsamis. Rawlings. Yates. &: Vitil (|2008l ) found that the fractional abundance of 



N2H + falls by an order of magnitude from cloud center to edge within a cloud undergoing inside-out 
collapse (see their Fig. 3). 

To account fully for the interplay of dynamics and chemistry, we would also need to follow 
simultaneously both a chemical reaction network and line transfer in our evolving cloud model. 
Needless to say, such an undertaking is beyond the scope of this study. However, we may assess in 
a general way the effect of abundance gradients. Cons ider the optically thick lines of HCN, w hose 



profiles should be especially sensitive. In the spirit of ITsamis. Rawlings. Yates. fe Vitil ([2008), we 



adopt a linear variation of the relative HCN abundance: 

X = X [1 + $ (M r /M - 1/2)] . (6) 
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Here, M r is the mass contained within any radius r and X is the nominal abundance we used 
previously. Within our canonical model, we vary the nondimensional slope (3 and calculate the 
emergent HCN profiles. 

The slope f$ cannot exceed 2 in absolute value, since X would then become negative at either 
the cloud center or edge. For smaller, but substantial, values of |/3|, the most significant effect 
is on the depth of the self-absorption dip. Figure 6 shows how the optically thickest, F (2 — 1), 
line profile changes as fi ranges from -0.5 through (uniform abundance), to +0.5. According to 
the standard gauge of the relative velocity shift, the effect is minor. For /3 = —0.5, <5Vhcn is 6 % 
greater in magnitude than the canonical value, while the shift is 6 % less for (3 = +0.5. 

Nevertheless, the degree of self-absorption diminishes markedly for positive /3, i.e., when HCN 
is concentrated in the more slowly moving, outer portion of the cloud. To restore the dip to its 
original magnitude, the contraction age for the cloud would need to increase slightly from the 
canonical 1.0 Myr to 1.2 Myr. We conclude that, while coupled chemical and dynamical models are 
preferable in principle, such extensive effort is unnecessary for beginning a quantitative assessment 
of starless core properties. 



Discussion 



This suite of radiative transfer calculations performed on our theoretical cloud model has 
yielded two important findings. First, we have been able to match the generic, observed proper- 
ties of the selected line profiles using a reasonable cloud mass, turbulent velocity dispersion, and 
evolutionary time. This success may indicate that the model has captured the essential physics 
of the contraction process itself. Note also that the accelerating character of the contraction in 
our theory explains naturally th e key observation that starless cores of lower density ex hibit less 
asymmetry in their line profiles (JGregersen &: Evansll2000i : IWilliams. Lee. fe Mversll2006l ). Second, 
the sensitivity of the emergent line profiles to the input parameters means that the model offers 
a practical route to gauge starless core properties that have previously been difficult to ascertain. 
For example, given the empirical mass and kinetic temperature of an object, the matching of line 
profiles would yield a first estimate of its contraction age. Obtaining ages for numerous starless 
cores in different environments would aid greatly in understanding the onset of star formation. 

Of course, the observations offer more detail than the model can explain, at least in its current, 
highly idealiz ed form. We cannot yet account for the rather surprising sp atial extent of asymmetric 



CS emission ( Myers. Evans. fe Ohashi 



2000 



Lee. Myers. fe Tafallall200ll ). Nor can we explain the 



relative intensities of the HCN hyperfine lines (ISohn et al.l 12004 ) . Explaining both features of 
the data may require alterations in the spatial abundance gradients, in the manner employed in 
Sections 4.2 and 4.4. Additionally, an improved model would consider more carefully how the 
interior, relatively quiescent gas joins onto surrounding, more turbulent, material. In fact, it has 
long been realized that the pattern of HCN hyperfine anomalies requires the presence of an envelope 
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of rela tively high excitation temperature ([Gonzales- Alfonso &: Cernichardll993l ). iGong fc Ostriker 
(J2009J) have recently simulated the birth of dense cores out of spherical, converging flows. Such 
calculations, when extended to three dimensions, may offer additional insight into the nature of 
the core-envelope transition. 

Finally, we must face the obvious fact that real starless cores are not spherical, as the model 
currently assumes. The shapes of these objects may be partially determine d by the anisotropic pres- 
sure of their turbulent environments, as found in recent simulations (e.g. lOffner. Klein. &: McKee 
2008). Another import ant factor setti ng the gross core properties is the force associated with the 
internal magnetic field (JLi et al.ll2004l ). It would be instructive to study, using the method of Pa- 
per I, the onset of collapse in m a rgina lly stable, magnetostatic equilibria, such as those found by 
Tomisaka. Ikeuchi. fc Nakamural (|1988l ). Such an investigation would add a significant degree of 
realism to our picture of protostellar collapse. 



We are grateful to Floris van der Tak, one of the authors of RATRAN, for guiding us in the 
proper use of the code. Jorge Pineda, who is independently analyzing the line emission from cloud 
models, also provided much-valued assistance, as did Melissa Enoch. Throughout the project, S. S. 
was partially supported by NSF grant AST-0908573. 
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Fig. 1. — Profile evolution of the F±,F = 0, 1 — > 1,2 hyperfine line within the J = 1 — > rota- 
tional transition of N2H + . The profiles, viewed through the cloud center, are for a cloud of 3 M & 
with k = 5.5 x 10~ 12 g cm -1 s~ 2 . From bottom to top, the times are 0, 0.5, and 1.0 Myr. 



18 



CD 

s_ 

CD 
Q. 

E 

CD 

I- 

co 

c 
c 

CD 

+-> 

c 

< 



1.2 


i 


1 ' 


1 ' 1 ' 




1.0 






CS(2-1) 

1 i 
1 i 


— 


0.8 






1 i 
1 i 

! * '* 
1 * /i 

1 1 l\ 

1 1 I 1 
1 I / 1 


— 








1 \ f 1 
1 * ' 1 
1 * ' 1 

1 * * 1 

1 * "* 




_ 


0.6 












0.4 










— 


0.2 




i 
i 
i 
i 






— 






i / 
/ / 
i / 
/ / 




- 


00 


" 


1 _^ 


i N^_ i 





■1.0 



-0.5 0.0 +0.5 

Velocity V r ( km s" 1 ) 



+ 1.0 



Fig. 2. — Profile of the J = 1 — > rotational transition of CS, viewed through the cloud center for 
the canonical cloud model. The dashed curve is the profile using a uniform CS relative abundance 
of X (CS) = 4 x 10~ 9 . For the s olid profile, the molecule was centrally depleted, following the 



prescription of iTafalla et al.l (J2002J). 
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Fig. 3. — Profiles of HCN for the three hyperfine lines within the J = l->0 rotational transition. 
All profiles are calculated using the canonical model, and are viewed through the cloud center. 



20 



0.0 



CO 



o 
o 

CD 

> 



-0.5 " 



■1.0 - 



■1.5 " 



-2.0 



» 

\\ 

\l 

\i 

V ' 

V /' 

V /' 

V /' 

V //' 

V / ' 

\v / ' 

OS. / ' 



o 



1 



17 



Radius r (10 cm) 








1 



17 



Radius r (10 cm) 



Fig. 4. — Theoretical velocity and density profiles. The left panel shows the contraction velocity 
as a function of radius. The solid curve is the canonical model (e = 0.1, t = 1 Myr), while the 
dashed curve has e = 0.2 and t = 0.25 Myr. The right panel displays the density profiles for the 
same two models. 
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Fig. 5. — The effect of changing e and t on the three hyperfine lines of HCN. All profiles are 
calculated using e = 0.2 and t = 0.25 Myr, and are viewed through the cloud center. Note how 
the asymmetry is decreased with respect to the profiles of Figure 3. 
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Fig. 6. — The effect of abundance gradients on the F (2 — 1) HCN hyperfine line. The cloud model 
is the canonical one, with the addition of a linear HCN abundance gradient. The slope /3, as defined 
in equation (6) is given in each panel. Notice how the self-absorption dip decreases for positive (5. 



